Anomaly at Ak in the angle-resolved photoemission spectrum of dirty superconductors 
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Elastic forward scattering can lead to an anomaly in the angle-resolved photoemission spectrum 
of the cuprate superconductors. Here we discuss how this anomaly can be used to provide a mea- 
surement of the superconducting gap Ak for k values away from the Fermi surface. 



In the search to understand the mechanism responsible 
for pairing in the cuprate high T c superconductors, a key 
signature is thought to be the momentum dependence of 
the gap. A gap which has the simple (cos k x — cos k y ) de- 
pendence throughout the Brillouin zone would imply that 
the pairing arose from a near neighbor Cu-Cu interaction 
such as a superexchange coupling. Additional k depen- 
dence, involving higher d x 2_ y 2 harmonics, would imply 
a more extended spatial interaction d], and some angle- 
resolved photoemission spectroscopy (ARPES) data on 
undcrdoped cuprates have indeed been used to suggest 
that the underdoped materials may have a longer range 
pairing^. Information on the gap function for values of 
k on the Fermi surface have been obtained from low tem- 
perature transport measurements 0, which probe the 
nodal regions, and from ARPES measurements^. While 
the peak at the quasiparticle energy E(k) in the ARPES 
energy distribution function (EDC) can in principle be 
used to determine A k for k values away from the Fermi 
surface, this requires assuming a band structure and ne- 
glecting self-energy shifts. In addition, the EDC can be 
broad and asymmetric because of interactions, making 
the determination of the peak position uncertain. Here 
we discuss an alternative possibility for determining the 
k dependence of the gap away from the Fermi surface 
for general values of k based upon a structure in A(h, uj) 
introduced by forward elastic scattering processes. 

One expects that in I^S^CaC^Os (BSCCO) out-of- 
plane impurities or disorder will lead to forward elastic 
scattering 0. A general discussion of the effect that such 
scattering can have on the ARPES spectrum was previ- 
ously given |(| 0| • Here we will focus on one aspect of the 
scattering which leads to an anomalous structural fea- 
ture in the ARPES spectral EDC at an energy equal to 
Ak- To illustrate the idea, we begin by considering just 
the effect of strong forward elastic scattering where the 
Nambu self-energies can be approximated by 



E (k,o;) = -ro(k) 



(i) 
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The conditions for this approximation to be valid were 
discussed in 0. Here r (k) is the normal state scat- 
tering rate due to out-of-plane impurities. Note, that 
expressions and possess square-root divergences 
at the gap edge u = ±Ajc. Using these self-energies, the 
one-electron spectral weight can be written as 



A(k,ui) = Im{-^ 
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where Z(k,ui) = 1 
following we will set 



zr (k) sgnw 



AT- 



In the 



£k = — 2i(cos k x + cos k y ) — At' cos k x cos k y — /i (5) 

with t' jt = —0.35 and \ijt = —1, giving the Fermi surface 
shown in Fig 1. We will also assume for the purpose of 
illustration that 



Ak = —^-{cos k x — cos k y ) 



(6) 



with An = 0.2t. 

Fig. 2(a) shows a plot of A(k,u>) for several values of 
the normal state forward scattering rate Fo(k) for k = 
(0.647T, 0) [the start of the A-cut shown in Fig. 1]. As 
expected, there is a broadened quasiparticle peak at lu = 
— \J ^ + A 2 .. However, in addition there is a square-root 
anomaly as to approaches — Ak- From eq. (@J) one finds 
that for k away from the Fermi surface with |ek| > Ak, 



A(k, 



-A k ) 
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(7) 



Thus, while the strength of the square-root anomaly 
caused by the forward scattering decreases as one goes 
deeper below the Fermi surface, the structure remains, 
and should be observable if the scattering rate To(k) is 
sufficiently large. 

In general, out-of-plane impurities will give rise to a 
scattering potential characterized by a finite range 
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and a more complete, self-consistent treatment of the self- 
energy due to impurity scattering is required. Here, we 
follow the approach taken in [(| in which a simple ex- 
ponential scattering potential Vbe _Kr was treated within 
a self-consistent Born approximation. In this case, the 
scattering strength depends upon both k and uj. For 
uj = — Ak, the scattering strength is suppressed when 
|ek| is greater than e K = vpK due to phase space restric- 
tions. Nevertheless, as shown in Fig. 2(b), the anomaly 
continues to occur at u — — Ak, although the structure 
can be broadened somewhat. Here, we have set k = 1, 
measured in inverse units of the Cu-Cu spacing, and cho- 
sen the strength of the out-of-plane impurity scattering 
potential and the impurity concentration to fix the nor- 
mal state scattering rate Tq equal to 2Ao, Ao, and Ao/2. 

In general, many-body interactions give rise to a gap 
function which depends upon both k and to. However, 
at low temperatures when to is equal to the real part of 
the gap at the gap edge i.e. Ak = Re[A(k, uj = Ak)], the 
broadening due to inelastic scattering vanishes as T 3 
and the imaginary part of the gap is determined by elastic 
impurity scattering only, i.e. the anomaly at Ak is stable 
against many-body interactions. Thus, out-of-plane for- 
ward impurity scattering gives rise to the anomaly and in- 
plane isotropic impurity scattering leads to only a small 
broadening for low energies near the Fermi level. In our 
units, a typical value for the broadening due to in-plane 
scattering is of order .02t, a typical value of the in-p lane 
scattering estimated from fits to STM studies It 
is irrelevant for the anomaly at the gap edge and other 
features at higher binding energies. 

In order to illustrate the type of behavior one is look- 
ing for, in Fig. 3 we have plotted A(k, uj) versus w for a 
set of k values taken along the A-cut shown in Fig. 1. 
Here, n = 1 and results including a self-consistent Born 
approximation treatment of the forward scattering, to- 
gether with an in-plane constant scattering rate of 0.02i, 
are shown for Tq = 2Ao, Ao, and Ao/2. In Fig. 3, the 
solid red curve follows the Ak anomaly while t he dashed 
red curve indicates the quasiparticle energy — ye^+ ^k- 



FIG. 2: a) Approximate spectral function A(k,ui) from Eq 
(4) versus to for k = (0.647T, 0) with To = 2Ao (dashed curve), 
Ao (solid curve) and Ao/2 (dotted curve), with Ao = 0.2t. b) 
shows the numerical result obtained from the self-consistent 
Born approximation as described in Ref. ||, with k — 1 
and the same scattering rates. Curves are offset by 0.5 with 
respect to each other and uj is given in units of t. 



The various A(k, uj) spectra are offset for the different 
k values which are taken at a momentum separation of 
Ak = 0.035 along the A-cut shown in Fig. 1. 

As seen, the strength of the anomaly depends 
upon the out-of-plane impurity scattering rate To. For 
To > Ao, there is a clear anomaly which occurs at 
lu = — Ak. Here, we have taken Ak = ^ (cos k x — cos k y ) 
for illustration. In general, the anomaly occurs at Ak = 
Re[A(k, oj = Ak)] and the locus of this anomaly allows 
one to determine the k-dependence of the gap for general 
values of k. 

Further results for A(\s.,uj) taken along the three mo- 
mentum cuts A, B, and C are illustrated in Fig. 4. Here 
k = 1 and Tq has been set to 2Ao. As before, the 
solid red curve gives the locus of — Ak, while the dashed 
red curve follows the broadened quasiparticle dispersion 
— \A k + ^k- ^ s seen a l° n § tne B cut, the strength of the 
gap anomaly decreases as Ak decreases. It is also sup- 
pressed when | ek | > e K due to phase space restrictions 
associated with the forward elastic scattering process. 

In Fig. [S] we show A(k, uj) at the M point k = (tt, 0) 
for three different values of k in order to illustrate the 
influence of the range of the scattering potential on the 
anomaly. A rough criterion for the appearance of the 
anomaly is that k should become smaller than 1, which 
means that the range of the scattering potential should 
become larger than one lattice spacing of the CuC>2 
planes. This is likely to be the case for out of plane 
impurities in the cuprates because of their poor screen- 
ing. 

In the foregoing analysis of the anomaly at Ak we have 
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FIG. 3: a-c A(\t,u>) versus to for k along the A cut shown 
in Figure 1 with a momentum separation Ak = 0.035. These 
plots illustrate how the anomaly varies for ft = 1 and different 
values of the the scattering strength: (a)To = 2Ao; (b) To — 
A ; (c) To = A /2, with A = 0.2t. The solid red curve 
indicates the Ak anomaly, while the dashed red curve follows 
the quasiparticle dispersion. Fermi surface crossing spectra 
are colored green. 



FIG. 4: a-c AQi, to) versus lo for k taken at intervals of 
5k = 0.0385, 0.12, 0.065 along the A, B, C cuts shown in Fig- 
ure 1, for Fo = 2Ao and ft = 1 . The solid red curve shows the 
locus of the anomaly which follows Ak = ■^ Q -(cosfe E — cosfcy). 
The dashed red curve follows the quasi-particle dispersion 
— v/ej + A£, and the green curve represents the EDC at the 
Fermi surface crossing point. 



restricted ourselves to the self-consistent Born approxi- 
mation. It would be interesting to study effects beyond 
Born approximation, for example by employing the T- 
matrix approximation as in a recent work by Rieck et 
al 01 > to see how this might affect the strength of the 
anomaly. The basic conclusion of our work, i.e. the ex- 
istence and position of the anomaly at lo = — Ak, will 
not be affected by the overall magnitude of the impu- 
rity potential, however (see, e.g. Eq. Q ). Wc leave the 
full calculation of the effect on ARPES of finite range 
scattcrers of arbitrary potential for a future study. 

In summary, the enhancement of the k-dependent den- 
sity of states at the gap edge lo = ±A k leads to an onset 
anomaly in A(k, lo) due to forward elastic impurity scat- 
tering. The locus of this anomaly provides a direct mea- 
sure of the k-dependent gap Ak at the gap edge. Strictly 
speaking, Ak = ReA(k, lo = Ak) where A(k, lo) is the 
complex k and w-dependent gap function. As noted, at 
low temperatures with co = Ak, the imaginary part of the 
gap arising from inelastic scattering vanishes as (T/Ao) 3 
and in-plane clastic scattering leads to only a small imag- 
inary contribution. Since out-of-plane forward elastic 
scattering gives rise to this anomaly, one would like to be 
able to control the number of surface defects, adjusting 
Tq to obtain an optimal measurement. One possibility 
may be to purposely add impurities to the surface. It 
is also possible that older BSCCO samples may contain 
sufficient out of plane disorder to see the effect discussed 
here with modern energy and momentum resolution. Al- 
ternatively, as the experiment proceeds and the surface 
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FIG. 5: A(k,u>) versus lo for k = (it, 0) with To = 2Ao and 
three different values of k = 2 (solid curve), ft = 1 (dotted 
curve) and ft = 0.5 (dashed curve), with Ao = 0.2t. 



gradually becomes contaminated, one may observe the 
development of this anomaly. 
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